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o 

H The geometry of a light wavefront evolving from a flat wavefront under the action of 

weak gravity field in the 3-space associated to a post-Newtonian relativistic spacetime, 
s ! is studied numerically by means of the ray tracing method. 
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1. Introduction 

' ' The curvature of initially plane light wavefronts by a gravity field is a purely general 

13^ relativistic effect that has no special relativistic analogue. In order to obtain an 

jj^ experimental measurement of the curvature of a light wavefront, Samuel recently 

W) proposed a method based on the relation between the differences of arrival time 

recorded at four points on the Earth and the volume of a parallelepiped determined 
by four points in the curved wavefront surface, see Ref. [T] In this work and in Ref. [2] 
with more detail, we study a discretized model of the wavefront surface by means 
of a regular triangulation for the study of the curvature(s) (mean and relative, see 



> 

o 
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Ref. 3) of this surface. 
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2. Light propagation in a weak gravitational field 
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Let us consider a spacetime (AA , g) corresponding to a weak gravitational field 
determined by a metric tensor given in a global coordinate system {(z,ct)} by 
9a/3 = Vap + h a fj, with rj a p = diag(l, 1, 1, — 1). Where the coordinate components 
of the metric perturbation h a p are: 

h ab = 2c- 2 K\\z\\~ 1 6 ab , h ai = -4c-*K\\z\\- 1 Z a , h 4i = 2c- 2 k\\z\\-\ (1) 

here k := GM represents the gravitational constant of the Sun, located at Z a (t), 
and c represents the vacuum light speed. The null geodesies, z(t) — (z(t),t\ satisfy 
the following equations, see Ref. [4] 

— (/l44,fc^? + h a k,l — \hu,a)Z % 

-(c-% fcli - ±c-*h jk , t )z>t k z a , 
= g afj z a zP. 

The second equation is the isotropy constraint satisfied by the null geodesies. 
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3. Local approximation of the wavefront 

Let be Sq a flat initial surface far from the Sun, formed by points (zi, z 2 , — £) (with 
C > 0) in an asymptotically Cartesian coordinate system {z}. For the discretization 
of iSo a triangulation is constructed in such a form that each vertex is represented 
by a complex number of the set: 

V := {z — ai + a 2 uj + a 3 uj 2 \ ai,a 2 ,a 3 £ A, uj := exp(27ri/3)}, (2) 

The initial triangulation by V induces a triangulation on the final wavefront St ■ The 
evolution of a photon z := z(0) £ V with velocity z :— (0,0, c) in phase space 
u = (z,z) may be written as a first order differential system u = F(u,t). This 
determines a flow, z(t) = (p t (zo, Zq), in the 3-dimensional curved quotient space of 
Ai by the global timelike vector field d t associated to the global coordinate system 
used in the post-Newtonian formalism. For each time t, the flow (p t determines a 
2-dimensional curved wavefront iS t . 

To compute the curvatures of the wavefront surface corresponding to the mesh V 
at each inner vertex, we consider a 1-ring formed by the six vertices closest. For each 
1-ring, one obtains on the mesh V the image under the flow ip t ■ In a neighbourhood 
of the image point the wavefront can be approximated by a least-squares fitting of 
the data obtained as the quadric: 

y 3 = f(y\y 2 ) ~ |ai(y x ) 2 + a 2 yV + h,{y 2 ) 2 . (3) 
using adapted normal coordinates {y 1 }. 



Numerical integrator 



We apply the ray tracing method, see Refs. 2[5 to a tubular region of light wavefront 



region supposing a gravitational model generated by a static Sun, considered as a 
point. The mean and relative total curvatures, defined in Refs. 2j3 are computed 



at each inner vertex of mesh on the light wavefront surface St in the vicinity of the 
Sun, by the implementation of the following pseudocode: 

Data: u* := «, z* n ), n = 1, . . . N 
for n = 1 . . . N do 

U n := Taylor(t,it*) 

y n := NormalCoordinates(;z n ) 

for i = . . . 6 do 
y rH := Ring(y n ) 

end 

(01,02,03) = LeastSquares(y„.) 
1ab{x u ) ■= Metric(ai, a 2 ,a 3 ,a; n ) 
B = SecondFundamentalForm(a; Jl ) 
(Ai, A2) = Diagonalize(j3) 
(K rc \,H) = Curvature(Ai, A2) 

end 
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Fig. 1. Wavcfront surface and relative curvature (gray scale) deformed by a spherical gravitational 
field (a different scale is used for the vertical axis). 



In Figure [T] the surface St at the time when the wavefront arrives at the Earth 
is shown using a gray-scale to represent the relative curvature (note we have used 
a different scale on the Oz 3 -axis) . One sees in this figure that the absolute value of 
the relative curvature defined on St increases as the distance between the photon 
and the Oz 3 -axis, where the Sun is located, decreases. 
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